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Abstract —We consider the delay minimization problem in 
an energy harvesting communication network with energy 
cooperation. In this network, nodes harvest energy from nature 
to sustain the power needed for data transmission, and may 
transfer a portion of their harvested energies to neighboring 
nodes through energy cooperation. For fixed data and energy 
routing topologies, we determine the optimum data rates, 
transmit powers and energy transfers, subject to flow and 
energy conservation constraints, in order to minimize the 
network delay. We start with a simplified problem where 
data flows are fixed and optimize energy management at each 
node for the case of a single energy harvest per node. This is 
tantamount to distributing each node’s available energy over 
its outgoing data links and energy transfers to neighboring 
nodes. For this case, with no energy cooperation, we show 
that each node should allocate more power to links with more 
noise and/or more data flow. In addition, when there is energy 
cooperation, our numerical results indicate that, energy is 
routed from nodes with lower data loads to nodes with higher 
data loads. We then extend this setting to the case of multiple 
energy harvests per node over time. In this case, we optimize 
each node’s energy management over its outgoing data links 
and its energy transfers to neighboring nodes, over multiple 
time slots. For this case, with no energy cooperation, we show 
that, for any given node, the sum of powers on the outgoing 
links over time is equal to the single-link optimal power over 
time. Finally, we consider the problem of joint flow control and 
energy management for the entire network. We determine the 
necessary conditions for joint optimality of a power control, 
energy transfer and routing policy. We provide an iterative 
algorithm that updates the data flows, energy flows and power 
distribution over outgoing data links sequentially. We show 
that this algorithm converges to a Pareto-optimal operating point. 


Index Terms —Energy cooperation, energy harvesting, wireless 
energy transfer, optimal routing, resource allocation 


I. Introduction 

We consider an energy harvesting communication network 
with energy cooperation as shown in Fig. [T] Each node 
harvests energy from nature and all nodes may share a portion 
of their harvested energies with neighboring nodes through 
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energy cooperation ID. We focus on the delay minimization 
problem for this network. The delay on each link depends on 
the information carrying capacity of the link, and in particular, 
it decreases monotonically with the capacity of the link for 
a fixed data flow through it; see e.g., m eqn. (5.30)]. The 
capacity, in turn, is a function of the power allocated to the 
link, and in particular, it is a monotonically increasing function 
of the power, for instance, through a logarithmic Shannon 
type capacity-power relationship; see e.g., 0 eqns. (9.60) and 
(9.62)]. In addition, the delay on a link is a monotonically 
increasing function of the data flow through it, for a fixed link 
capacity l2] eqn. (5.30)]. 

In this paper, we consider the joint data routing and capacity 
assignment problem for this setting under fixed data and en¬ 
ergy routing topologies El Section 5.4.2]. Our work is related 
to and builds upon classical and recent works on data routing 
and capacity assignment in communication networks m, a- 
na, and recent works on energy harvesting communications 
na-ini and energy cooperation ID, na-Gi in wireless 
networks. In our previous work ID, Gi, we studied the 
optimal energy management problem for several basic multi¬ 
user network structures with energy harvesting transmitters 
and one-way wireless energy transfer. Inspired by joint routing 
and resource allocation problems in the classical works such 
as a-Q, Qoi, ca, in our current work, we study joint 
routing of energy and data in a general multi-user scenario 
with data and energy transfer. We specialize in the objective 
of minimizing the total delay in the system. To the best of 
our knowledge, this problem has not been addressed in the 
context of energy harvesting wireless networks with energy 
cooperation. Among previous works, the approach that is 
most related to ours is that in reference ll26l . which studies 
networkwide optimization of energy and information flows 
in communication networks with simultaneous energy and 
information transfer. We also note the references ED, m 
for related joint data routing and energy transfer schemes 
in networks with special energy transfer capabilities and no 
energy harvesting. Finally, we refer the reader to El-Ol 
for a related line of research about resource allocation in base 
stations powered by renewable energy and energy cooperation. 

We divide our development in this paper into three parts. 
In the first part, we assume that the data flows through the 
links are fixed, and each node harvests energy only once. In 
this setting, we determine the optimum energies allocated to 
outgoing data links of the nodes and the optimum amounts of 
energies transferred between the nodes. In the second part, we 
extend this setting to the case of multiple energy harvests for 
each node. In the last part, we optimize both data flows on 
the links and energy management at the nodes. We determine 
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Fig. 1. System model. 

the jointly optimal routing of data and energy in the network 
as well as distribution of power over the outgoing data links 
at each node. 

In the first part of the paper, in Section [Till we focus on 
the optimal energy management problem at the nodes with 
a single energy harvest at each node. First, we consider the 
case without energy cooperation. We show that this problem 
can be decomposed into individual problems, each one to be 
solved for a single node. We show that more power should 
be allocated to links with more noise and/or more data flow, 
resembling channel inversion type of power control ll^ . Next, 
we consider the case with energy cooperation, where nodes 
transfer a portion of their own energies to neighboring nodes. 
In this case, we have the joint problem of energy routing 
among the network nodes and energy allocation among the 
outgoing data links at each node. For this problem, we develop 
an iterative algorithm that visits all energy links sufficiently 
many times and decreases the network delay monotonically. 
We numerically observe that energy flows from nodes with 
lightly loaded data links to nodes with heavily loaded data 
links. 

In the second part of the paper, in Section HVl we extend our 
setting to the case of multiple energy harvests at each node, 
by allowing time-varying energy harvesting rates over large 
time frames. We incorporate the time variation in the energy 
harvests and solve for the optimal energy management at each 
node and energy routing among the nodes. First, we focus on 
the case without energy cooperation. We show that the sum 
powers on the outgoing data links of a node over time slots is 
equal to the single-link optimal transmit power of that node 
over time and can be found using na-iiBi. When the optimal 
sum powers are known, we show that the problem reduces 
to a problem with a single energy arrival and can be solved 
using our method. Next, we focus on the case with energy 
cooperation. We show that this problem can be mapped to the 
original problem with no energy cooperation by constructing 
an equivalent directed graph. 

In the last part of the paper, in Section |Vl we consider 
the problem of determining the jointly optimal data and 
energy flows in the network and the power distribution over 


the outgoing data links at all nodes. We determine a set 
of necessary conditions for the joint optimality of a power 
control, energy transfer and data routing policy. We then 
develop an iterative algorithm that updates the data flows, 
energy flows and distribution of power over the outgoing data 
links at each node in a sequential manner. We show that this 
algorithm converges to a Pareto-optimal operating point. 

II. Network Flow and Energy Model 

We use directed graphs to represent the network topology, 
and data and energy flows through the network. All nodes 
are energy harvesting, and are equipped with separate wireless 
energy transfer units. Information and energy transfer channels 
are orthogonal to each other. 

A. Network Data Topology 

We represent the data topology of the network by a directed 
graph. In this model, a collection of nodes, labeled n = 
1,...,/V, can send and receive data across communication 
links.In particular, a node can be either a source node, a 
destination node or a relay node. A data communication link 
is represented as an ordered pair (/, j) of distinct nodes. The 
presence of a link (i,j) means that the network is able to 
send data from the start node i to the end node j. We label 
the data links as / = 1,..., T. The network data topology can 
be represented by an x T matrix. A, in which every entry 
Ani is associated with node n and link I via 

{ 1, if n is the start node of data link I 
— 1, if n is the end node of data link I (1) 
0, otherwise 

We define Od{n) as the set of outgoing data links from node 
n, and Irf(n) as the set of incoming data links to node n. We 
define TV-dimensional vector s whose nth entry s„ denotes the 
non-negative amount of exogenous data flow injected into the 
network at node n. On each data link I, we let ti denote the 
amount of flow and we call the L-dimensional vector t the 
flow vector. At each node n, the flow conservation implies: 

ti- y] ti= Sn, Vn (2) 

lGOd{n.) lGXd{n) 

The flow conservation law over all the network can be com¬ 
pactly written as: 

At = s (3) 

We define c/ as the information carrying capacity of link 1. 
Then, we require ti < ci, VI. 

B. Network Energy Topology 

All nodes are equipped with energy harvesting units. In this 
section, we describe the energy model for the case of a single 
energy harvest per node. We present the extension to the case 
of multiple energy harvests in Section HV] Here, each node n 
harvests energy in the amount of We use TV-dimensional 
vector E to denote the energy arrival vector for the system. In 
the energy cooperation setting, there are energy links similar 
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to data links. An energy link is represented as an ordered pair 
{i,j) of distinct nodes where the presence of an energy link 
means that it is possible to send energy from the start node 
to the end node. Energy links are labeled as q = 1,... ,Q. 
Energy transfer efficiency on each energy link is denoted with 
0 < ttq < 1 which means that when S amount of energy is 
transferred on link q from node i to node j, node j receives 
aq6 amount of energy. We assume that the directionality and 
the position of energy transfer links are fixed whereas the 
amount of energy transferred on these links are unknown. The 
network energy topology can be represented by an x Q 
matrix, B, in which every entry Bnq is associated with node 
n and energy link q via 

{ 1, if n is the start node of energy link q 
—aq, if n is the end node of energy link q (4) 
0, otherwise 

On each energy link q, we let yq be the amount of energy 
transferred. We call the E-dimensional vector y the energy 
flow vector. We denote by Oe{n) and Xe(n), respectively, the 
sets of outgoing and incoming energy links at node n. 


C. Communication Model and Delay Assumptions 

Eollowing the M/M/1 queueing model in ||2l, we represent 
the delay on data link / as; 


Di 


ti 

Cl - ti 


(5) 


where // is the flow and cj is the information carrying 
capacity of link I, with ti < ci, VL This delay expression 
is a good approximation for systems with Poisson arrivals 
at the entry points, exponential packet lengths and moderate- 
to-heavy traffic loads m. In view of energy scarcity in the 
network, moderate-to-heavy traffic load assumption generally 
holds. The packet arrival and packet length assumptions are 
made for convenience of analysis. Moreover, we assume that 
the slot length is sufficiently large to enable convergence to 
stationary distributions. In particular, we assume that the slot 
length is sufficiently longer than the average delay yielded by 
the M/M/1 approximation. Each node n, on the transmitting 
edge of data link I, with channel noise ai, enables a capacity 
Cl by expanding power pi. These quantities are related by the 
Shannon formula (|3] eqn. (9.60)] as; 




( 6 ) 


where all logs in this paper are with respect to base e. At each 
node n, the total power expanded on data and energy links are 
constrained by the available energy, i.e.. 


Pi+ X] yq^^ri+ y] Vn (7) 

l€Od{n) q^Ocin) qGl^{n) 


Using i-dimensional vector p = {pi,... ,pl) and F = A+ 
where = max{A„j,0}, the energy availability con¬ 

straints can be compactly written as; 


Fp + By < E 


( 8 ) 


We note that we use power and energy interchangeably in ([8]l 
and in the rest of the paper by assuming slot lengths of 1 unit. 


III. Capacity Assignment Problem for Single Time 

Slot 


In this section, we consider the capacity assignment problem 
for the case of a single energy harvest per node. We assume 
that the flow assignments, ti, on all links are fixed and are 
serviceable by the harvested energies and energy transfers. The 
total delay in the network is; 


^ = E 

I 


ti 

Cl - tl 


(9) 


The capacity assignment problem, with the goal of minimizing 
the total delay in the network is; 


min 

Cl,Pl,Vq 

s.t. 


E' 


Cl - tl 

I 

Fp + By < E 

tl < Cl, V/ 


( 10 ) 


By using the capacities ci in (|6]l, we write the problem in 
terms of the link powers pi and energy transfers pq only as; 



s.t. Fp + By < E 

Pi>(Ti{e^^‘-l), V/ (11) 

We solve the problem in (fTTI) in the rest of this section. We 
first identify some structural properties of the optimal solution 
in the next sub-section. The following analysis relies on the 
standing assumption that this problem has at least one feasible 
solution. To see if this problem is feasible, one can replace 
the objective function of (fTTI) with a constant and solve a 
feasibility problem, which turns out to be a linear program. 


A. Properties of the Optimal Solution 

First, we note that the objective function can be writ¬ 
ten in the form Myi^i)) where and 

g{xi) = ^ log (1 + Xi). Since / is convex and non-increasing 
and g is concave, the resulting composition function is convex 
ll3^ . The constraint set is affine. Therefore, (dU) is a convex 
optimization problem. The Lagrangian function is: 



E Pi 

ieOd(") 


+ E 2^9- 

qGOe{n) 

- E h - cn (e^*‘ - l)] - E ^9^9 

I q 


^ ^ ^qVq 
q€leM 


where A„ and jdi are Lagrange multipliers corresponding to the 
energy constraints of the nodes in ([8]) and the feasibility con¬ 
straints tl < Cl, respectively. The KKT optimality conditions 
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are: 


h'i{,pi) + Xnii)-l3i=G, VI (13) 

^m{q) ^q^k{q) ^q 0 ; (14) 

where hi{pi) = h (^5 log , n{l) is the be¬ 

ginning node of data link I, m{q) and k{q) are the beginning 
and end nodes of energy link q, respectively. The additional 
complementary slackness conditions are: 

A« Pi + X] Vq- X “ 92/9 = 0, Vn 

\lGOd{n) q&Oe{n) q^I^{n) J 

(15) 

Pi [pi - (Ti - 1 )] = 0 , VI ( 16 ) 

0qyq = 0, Vq (17) 

We now identify some properties of the optimal power 

allocation in the following three lemmas. 

Lemma 1 If the problem in ( 1771 ) is feasible, then Pi = 0, VI. 

Proof: If the problem in (fTTT i is feasible, its objective function 
must be bounded. Equality in the second set of constraints in 
(fTTIi for any I implies that the objective function is unbounded. 
Therefore, we must have strict inequality in those constraints 
for all I, and from (fThl l. we conclude that Pi = 0,VI. ■ 

Lemma 2 At every node n, the optimal power allocation 
amongst outgoing data links satisfies 

K{Pi) = h'm{Pm), Vfm&Od{n) (18) 

Proof: From (fOl l and Lemma [T] we have, 

Kipi) = -K{i), V/ (19) 

For outgoing data links I and m that belong to the same node 
n, 

KiPl) =-Xn= h'^{pm) (20) 


B. Solution for the Case of No Energy Transfer 

In the case of no energy transfer, we have yg = 0, Vq, and 
the problem becomes only in terms of pi as stated below: 


min 

Pi 


X 


_ tj. _ 

1 log {1 + ^)- ti 


s.t. X 

(eOd(n) 

Pi > cri (e^‘‘ - 1 ) , VI (25) 

This problem can be decomposed into N sub-problems as: 


min 

PI 

s.t. 


X X 

n leOd(n) 



^ Pi < En, Vn 

l&Od{n) 

Pi > cri(e'^*‘ - 1), VI 


(26) 


Since the constraint set depends only on the powers of node 
n, there is no interaction between the nodes. Every node will 
independently solve the following optimization problem: 


min 

Pi 

s.t. 


X 

leOdin) 



X Pi < En 

leOdin) 


Pi > cri - 1 ) , VI e Od{n) (27) 

The feasibility of ( |27l i requires En ^ J2lFOdln) 1) 

which we assume holds. Similar to (HU), ( | 2 U is a convex 
optimization problem with the KKT optimality conditions: 


h'i{pi) + X = 0, Vl€Od(n) 
with the complementary slackness condition: 


A I X Pi - En] = 0 

yleOd{n) 


(28) 


(29) 


which gives the desired result. ■ 

Lemma 3 If some energy is transferred through energy link 
q across nodes (i,j), then, 

K(pi) = Ctqh'miPni), VI & Od{i) , Vm € Od{j) (21) 

Proof: If some energy is transferred through energy link q, 
then yg > 0, and from (fTTl i. 9g = 0. From (fl4l i. we have, 

Xi = UgXj (22) 

Writing (fT3l l for nodes i and j, we have, 

h'liPi) = -A„ VI e Odd) (23) 

h'm{Pm) = -Xj, Vm e Od{j) (24) 

and the result follows from combining (|22] |. (l23b and (l24l i. ■ 

In the following two sub-sections, we separately solve the 
problem for the cases of no energy transfer and with energy 
transfer. 


The Lagrange multipliers for the second set of constraints in 
dZTl ) are not included, because similar to Lemma [U they will 
always be satisfied with strict inequality. From ( 1281 ), we have 


-A = h'fipi) 


-ti 

2ai 


ilog 

7 + -) 

2 





(30) 

(31) 


After some algebraic manipulations shown in Appendix lAl we 
have 


(32) 


where = \/VS- and W{-) is the Lambert W function 
defined as the inverse of the function w —>■ we^ EU- Next, we 
prove some monotonicity properties for the optimal solution, 
as a function of the qualities of the channels and the amounts 
of data flows through the channels. 


Lemma 4 For fixed ti, pi is monotone increasing in ai. 
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Proof: By differentiating (l32l i and using the following prop¬ 
erty 1^ 


dW{x) W{x) 

dx x{l + W{x)) 

it can be verified as shown in Appendix |B] that 


dpi 

dai 


„2t, 


l + Wizi) 


- 1 > 0 


(33) 


(34) 


where the inequality follows from > 1, Vf; > 0, and 
Ypj- > 1, Vz > 0, proving the lemma. ■ 

This lemma shows that, for fixed data flows, more power 
should be allocated to channels with more noise power, similar 
to channel inversion power control ll^ . 


Lemma 5 For fixed ai, pi is monotone increasing in ti. 


Proof: By differentiating (l32l i. it can be verified as shown in 
Appendix |C] that 

dpi ai{W{zi) + 2ti)e^^^(^d+U) 

dti ti{l + W{zi)) 

proving the lemma. ■ 

This lemma shows that, for fixed channel qualities (i.e., 
fixed noise powers), more power should be allocated to links 
with more data flow. 

Finally, we solve as follows: From the total energy 
constraint, we have J^iPii^*) — We perform a one 
dimensional search on A to find A* that satisfies J^iPii^*) ~ 
En, where pi{X*) is given in ( |32] |. Once A* is obtained, the 
optimal power allocations are found from (l32l i. 


C. Solution for the Case with Energy Transfer 

Now, we consider the case with energy transfer, i.e., yq > 0 
for some q. Assume that some energy pg > 0 is transferred 
from node i to node j on energy link q. Writing (l32l i for the 
outgoing data links of node i and node j, we have, 

Pi(Az) =a; -l) , ^l£OS) (36) 

V/eOd(j) (37) 

where za = and z^i = From (EDi, we 

have Xi = aqXj. The energy causality constraints on node i 
and node j are: 

^ pi{X*)=E,-yg (38) 

Fi(A*) = Ej + aqPq (39) 

ieOdO) 

Equations (l22l i. (l38l l and (|39] | imply 

o-q X! Pii'^q^*j)+ X] Pi{X*) = oiqEi + Ej (40) 

l&Od{i) l&OdU) 


Algorithm 1 Algorithm to solve capacity assignment problem 
for single time slot 

Initialize > No energy transfer 

1 : for I = 1 : W do 

2 : Find Ai such that Y.i^oi(i) Pi^^i) = Pi •ESll 

3: end for 

Main Algorithm 

4: for g = 1 : Q do t> All energy links 

5: Set (i,j) ^ (origin,destination) of energy link q 

6 : if Xi < aqXj then i> Perform energy transfer 

Find A* such that 

3 ” ’^leOdU) ~ i^q^i + Ej 

Set Tap^ = Ei- J^ieOdii) Pi^^q^]) ^ Update tap level 

0 Update battery levels 

Set Ei = X]jGOd(i)-Uj = J2ieOd{j)Pi^^j^ 

7: else if Xi > aqXj then i> Recall some energy 

8 : while Tapg > 0,Ai > aqXj,Ej > 0 do 

0 Recall e energy 

Set Ei = Ei + e, Ej = Ej - age, Tap^ = Tap^ - e 
Find Xi , Xj such that 

9: end while 

10 : end if 

11 : end for 


often. Since we do not know which energy links will be 
active in the optimal solution, we may need to call back any 
transferred energy in the previous iterations. To perform this, 
we keep track of transferred energy over each energy link by 
means of meters as in IT]. Initially, we start from the no energy 
transfer solution and compute A„ for every node n as described 
in the previous section. At every iteration, we open only one 
energy link g at a time, and whenever energy flows through 
link g, (l40l i must be satisfied with Ei and Ej in (l40l i replaced 
with the battery levels of nodes i and j at the current iteration. 
In particular, if Xi < aqXj, we search for A* that satisfies (l40l) . 
If no solution to (l40l) can be found, this means Ai > aqXj, 
and then previously transferred energy must be called back 
to the extent possible according to the meter readings. The 
algorithmic description is given below as Algorithm 1. From 
the strict convexity of the objective function, we note that each 
iteration decreases the objective function as described similarly 
in m Section V.A]. Our algorithm converges since bounded 
real monotone sequences always converge, and the limit point 
is a local minimum because, the iterations can only stop when 
Xi = aqXj for the energy links where pq > 0 which are the 
KKT optimality conditions from (l22l i. This local minimum is 
also the unique global minimum due to the convexity of the 
problem. 

IV. Capacity Assignment Problem eor Multiple 
Time Slots 


which can be solved by a one-dimensional search on A*. 

We solve (fTTT i by iteratively allowing energy to flow through 
a single link at a time provided all links are visited infinitely 


In this section, we consider the capacity assignment problem 
for the scenario where the energy arrival rates to the nodes can 
change over time. We assume that the time is slotted and there 
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are a total of T equal-length slots. In slots i = 1,..., T, each 
node n harvests energy with amounts Eni,En 2 , • ■ ■, EnT, and 
the arriving energies can be saved in a battery for use in 
future time slots. The subscript i denotes the time slot, and 
the quantities tii,cii,pu,aii and y^i denote the flow, capacity, 
power, noise power, and energy transfer in slot i. We assume 
that the flow allocation and channel noises do not change 
over time, i.e., tu = ti and an = ai,yi,yi. We further 
assume that the slots are long enough so that the M/M/1 
approximation is valid at every slot i. In particular, slot length 
is sufficiently larger than the average delay resulting from the 
M/M/1 approximation. Then, the average delay on link I at 
time slot i is given as. 


Dii 


ti 

tl 


(41) 


where cu = ilog^l + ^j. As the energy that has not 
arrived yet cannot be used for data transmission or energy 
transfer, the power policies of the nodes are constrained by 
causality of energy in time. These constraints are written as: 


E 


E E 

l&Od(n) q^Ocin) 

k / 



< 


'y ^ + y ^ CXqUqi j j 

^ —1 \ (jGXe(n) / 


Vn, V/c 


(42) 


The capacity assignment problem with fixed link flows to 
minimize the total delay over all links and all time slots can 
be formulated as: 


min 

Pli,Vqi 


s.t. 




EEt 

I 2 

E f E E 2^9*) 

j=l \ieC>d(ra) q^Oe{n) / 

k / \ 

— y ' ( P^ni + y ^ Olqyqi 

i=l \ ijGle(n) 

PH > crz(e^*' - 1), V/, Vi 


Vn, Vfc 


, V / c/, 


(43) 


The problem in (l43t is convex and the Lagrangian function 
can be written as: 


^=EE hiipii) + EE A- 


2=1 I 


n 


E Up- 

2=1 \l^Od{n) 


“1“ ^ ^ Vqi Ejii ^ ^ (^qUqi 

q^Oe{n) gGXe(n) 

T 

-EE (^qiUqi 


(44) 


q i—1 


where 


. The Lagrange 


hiipii) = ti |log(l + ^) -ti 
multipliers for the second set of constraints for ( |4^ are not 
included here because similar to before, they will always be 
satisfied with strict inequality. The KKT optimality conditions 


are: 

T 

K{ph) + X! ^nii)k = 0, V/, Vi (45) 

k=i 

T T 

^ ^ ^7n{q)k ^q ^ ^ ^r{q)k ^qi ~ 0 ; (^ 6 ) 

k—i k—i 

where n(/) is the beginning node of data link /, m{q) and 
r{q) are the beginning and end nodes of energy link q. The 
additional complementary slackness conditions as: 



E E Pli T y ^ yqi Ejii 

i=l \lGOd{n.) q&Oe{n) 


/ , CtqUqi 
qGle{n) / 


= 0 , 


Vn, V/c 


^qiVqi — O 5 V(^, Vi 


(47) 

(48) 


Now, we extend Lemmas |2] and [3 to the case of multiple 
energy arrivals over time. 


Lemma 6 Af every node n, the optimal power allocation 
amongst outgoing data links satisfies 

h'liPh) = h'^{pmi), Vl,m G Od{n),Vi (49) 


Proof: From (l45l i. we have, 

T 

h'liPli) =-'^K{i)k (50) 

k—i 

For outgoing data links I and m that belong to the same node 
n, 

T 

^liPli) ~ ^ ^ ^nk — ( 51 ) 

k—i 

from which the result follows. ■ 


Lemma 7 If some energy is transferred through energy link 
q across nodes (a, b) at time slot i, 

h'fipii) = aqh’^{prm), V/ S Od{a), Vm e Od{b) (52) 

Proof: If some energy is transferred through energy link q at 
time slot i, then yqi > 0, and from ( l48l l, 6qi = 0. From ( l46b . 
we have, 

T T 

y ^ ^ak — Olq y ^ Xbk (53) 

k—i k—i 

Then, we have, 

T T 

^liPli) ~ ^ ^ ^ak — Ckq ^ ^ ^bk 

k—i k—i 

= aqh'^{pmi), V/ e Od{a), Vm £ Od{b) (54) 

where the first equality follows from writing ( |45] | for node a, 
the second equality follows from ( |53] ), and the third equality 
follows from writing (|45]) for node b. ■ 
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In the following two sub-sections, we separately solve the 
problem for the cases of no energy transfer and with energy 
transfer. 


A. Solution for the Case of No Energy Transfer 


In this case, we have yqi = 0,Vi, Vg. The problem becomes 
only in terms of pu as follows: 

min 

Pii 

YY ' 

i=l 1 


s.t. 

k k 

E E Pli<'^Erii, Vn,\/k 

i=l leOdin.) i=l 



Pli > - 1), V(, Vi 

(55) 

The problem can be decomposed into N sub-problems 

as: 

min 

PH 

YY Y ' 

^ n ze^„)ilog(l + ff)-f, 


s.t. 

k k 

E E pii<^Eni, yn,yk 

i=l teOdin) i=l 



Pli > criY*'- yi,yi 

(56) 


Since the constraint set depends only on the powers of node 
n, there is no interaction between the nodes. Every node will 
independently solve the following optimization problem: 


i=l l&Od{n) 

tl 

ilog(l + ff)-f, 

k 

k 

E E 

Pli — ^ni: V/c 

*=1 leOdin) 

i=i 

Pli > <JiY 

'-1), \/lGOdin), 


(57) 


Solving ( ISTT i entails finding the optimal energy management 
policy for each link I, over all time slots i. We define bu = 
Pu-cri{e'^^‘ - 1) and G„i = E^i - \Od{n)\ai{e'^*‘ - 1). Then, 
(Ell becomes: 


min 

hi 


s.t. 


T 


E E 

iGOd(n) 



k 


k 


E E 

i=l leOd{n) i=l 


V/c 


bu >0, W G Odin), 'ii 


(58) 


For feasibility of ( fSSl l we need G„i > 0 which we assume 
holds. Now, we state an important property of the optimal 
policy which is proved in Appendix iDl 


Lemma 8 The optimal total power allocated to outgoing data 
links at each slot i, 'Yhi^Odin) same as the single-link 

optimal transmit power with energy arrivals G„i. 

From Lemma [8] we have that the sum powers in outgoing 
data links are given by the single-link optimal transmit powers 


which can be found by the geometric method in lfT3l or by the 
directional water-filling method in US). Once the sum powers 
are obtained, individual link powers are found by solving xisf 
which is defined in ( |92] ) in Appendix iDl The problem in x{si) 
includes a single energy harvest and is in the form of dZTI) . 
therefore, we use the method proposed in Section llll-BI to find 
the individual link powers. 


B. Solution for the Case with Energy Transfer 

From (l45l) and some algebraic manipulations we have 

Pu = m - l) (59) 

where zu = , ‘ —r— and W{-) is the Lambert W 
function. The Lagrangian structure of this problem is more 
complicated compared to the previous case since the power 
allocation at time i depends on {^n{i)k]"k=i- Therefore, here, 
we offer an alternative solution. 

In the scenario described above, the nodes have the capa¬ 
bility to save their energies to use in future slots. We note that 
saving energy for use in future slots is equivalent to transfer¬ 
ring energy to future slots with energy transfer efficiency of 
a = 1. In light of this observation, an equivalent representation 
of (l43l l can be obtained by modifying the network graph where 
each time slot is treated as a new node with a single energy 
arrival and saving energy for future slots is represented by 
energy transfer links of efficiency 1. The modification to the 
network graph is performed in the following way. First, we 
make T replicas of the network graph including all the nodes 
and the existing data and energy transfer links. Each replica 
will denote the network at one time slot. We let each replica 
node receive one energy harvest which amounts to the energy 
harvested by that node in that time slot. We keep the existing 
energy and data links but we add new energy links between 
different replicas of the same node. For every node n, we add 
energy links of efficiency 1 between replicas k and where 
k = 1,... ,T—1. Relabeling the nodes, we obtain a new graph 
where all nodes have one energy harvest. Essentially, we have 
reduced this problem to the case in Section IIII-CI and we use 
the solution provided in that section. 

We finally remark that our framework can easily be ex¬ 
tended to address variations in channel fading coefficients and 
energy transfer efficiencies by allowing the noise powers ai 
and energy transfer efficiencies ai to vary from slot to slot, 
i.e., defining cu = ^ log (1 + replacing ai with au- 


V. Joint Capacity and Flow Optimization 


In this section, we consider the joint optimization of capac¬ 
ity and flow assignments, in contrast to capacity assignment 
only with fixed flows, as considered in the previous sections. 
We focus on the case with a single energy harvest per node 
as in Section |III] The delay minimization problem with joint 
capacity and flow allocation can be formulated as: 


min 

PUVqM 


S.t. 


E 


ti 


Fp -f By < E 













(60) 


(67) 

( 68 ) 

(69) 


Pi > o-z(e^*‘- 1), V; 

At = s (60) 

where we optimize not only the powers pi and energy transfers 
Dq, but also the data flows ti. In ( l60b . the first set of constraints 
are the energy constraints, the second set of constraints are 
the capacity constraints on individual links, and the last set of 
constraints are the flow conservation constraints at all nodes. 

We assume that the exogenous arrivals s is serviceable by 
the energy harvests and energy transfers. This means that 
problem dflOb has a bounded solution and furthermore no data 
link is operating at the capacity, i.e., the capacity constraints 
are never satisfied with equality unless ti = pi = 0. We 
solve the problem in dhOl l in the remainder of this section. 
Here, the constraint set is convex, however, the objective 
function is not jointly convex in pi and ti 12, therefore, ( l60t 
is not a convex optimization problem. We study the necessary 
optimality conditions by writing the Lagrangian function as 
follows; 


"-St 


ti 






+ X! Vq-^n- X! “92/9 

q£Oe{n) q£le{n) 


E 

ieOd(n) 

- E 


- 1 )] +X 1 

n 

- E ^ 92/9 - E 7*^* 

q I 

The KKT optimality conditions areQ 


E E 

l&Odiji) leZd{n) 


-tl 

2ai 


-A = o, V/ 
1 


1 + ^ 

cri 


+ A. 


n{l) 






1 -2 


- tl 


+ Vn(i) - Vm{i) -11 + 2Piaie ' = 0, VZ 
Afc(q) - aqX^tq) - 6»q = 0, V(? 


Vr,. \ ^ tl- ^ j = 0, Vn 

lieC>d(n) leXd{n) 


PqVq = lltl = 0, Vg, yi 

I3i[pi - (Ji{e'^*‘- 1)] = 0, VI 
An, A, Oq, 11 > 0, VI, Vq, Vn 

We note that i/n < 0 is allowed since the Lagrange multiplier 
n corresponds to an equality constraint. Lemma |9] proved in 
Appendix [E] states the necessary optimality conditions. 

Lemma 9 For a feasible set of flow variables {ti}fli, 
transmission power allocations {pi}fli and energy transfers 
{yq}^=i to be the solution to the problem in ( I6QD . the following 
conditions are necessary. 

1) For every node n, there exists a constant An > 0 such that 

1-2 / X -1 

1 + ^ 

VI G Od{n) 


tl 

2ai 


1 


log ( 1 + ^ ) - fz 


S An, 


(70) 


and with equality if pi > 0. 

2 ) For every node n, there exists a constant bn > 0 such that 

-2 


E ' + 




EL 

O'; 


Vd=l,...,D 


- tl 


(71) 


(61) 


where Tn,d is a data path that starts from node n and ends 
at destination node d and for which pi > 0,VI G LFn,d- The 
condition in ([73 is valid for all such data paths that start 
from node n and end at any destination node. 

3) For all energy transfer links q, and VI G Od{n),Vk G 
Oditn) such that pi > 0 and pk > 0 where n and m are the 
origin and destination nodes of energy transfer link q 


tl 

2ai 


1 


log ( 1 + ^ ) - fz 


cri 


(62) 


(63) 

(64) 


CXq 


'2ak 


1 


■ log 1 + 




1 + ^ 
<Ji 

1 -2 
tk 


> 

Ctk 


(72) 


where m is satisfied with equality if Pq > 0. 


where nil) and m{l) are the source and destination nodes of 
data link I, k{q) and z{q) are the source and destination nodes 
of energy link q, respectively. The complementary slackness 
conditions are: 

An E yq~ ^ aqPq = 0, Vn 

\t£Od{n) qGOein) qele{n) J 

(65) 


From Lemma |9] the structure of the optimal solution is 
as follows: We define hi{pi,ti) as the objective function of 

the problem in (|60ll, hfipuf) = f A log (^1 + - fz 

We see from dTOl i that nodes should allocate more power on 


links where the quantity 


dhi 

dpi 


is large and less power on 


( 66 ) 


^With the objective function of there is an uncertainty when ti = 
Pi = 0. Nonetheless, we argue as in (2 page 441] that the objective function 
of <60t is differentiable over the set of all pi with ^ log > ti and 

= 0, = 0 and = 0 are necessary conditions for optimality. 


links where this quantity is small. Similarly, from ( ItTI ). we see 
that less flow should be allocated on paths where the quantity 
^ is large and more flow on paths where this quantity 

is small. Finally, ( |72] ) tells us the necessary conditions for 
energy transfer. We describe our solution to the problem in 
dflOl) in the next section. 

A. Algorithmic Solution for the Joint Capacity and Flow 
Optimization Problem 

In this section, we propose an iterative algorithm. There are 
three steps to each iteration as summarized below. We start 
from a feasible point (to,po). 
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1) Energy Management Step: We fix a stepsize > 0- 
Each node computes ^ for their own outgoing data 
links where pi > 0. Every node performs the following 
iteration; 


= 


\Pi if ^ = arg max,gc)^(„) 

\Pi if ^ = arg minigci_j(„) 


Oh, 


dpi 
dhi 


dpi 


VPl 


Otherwise 


(73) 


where k denotes the iteration number, and the derivatives 
are computed at the current iteration, i.e., for 
2) Data Routing Step: We fix a stepsize > 0. Each node 
n computes ^ ^ for the data paths originating 

l&tFn,d 

from source node n and ending at any destination. 
Assume the path T* maximizes ^ ^ and the 

n,d 

path Q* minimizes ^ ^ for each n. Every node 

performs the following iteration: 

rff-6, tfi&K 

+ if^eg: (74) 

[ tf , otherwise 


3) Energy Routing Step: This step is the same as described 
in Section IIII-CI Specifically, every node goes through 
its energy transfer links and makes the comparison 
where m is the receiving node of 


dhi 

dpi 




dh-n 


energy link q. If 


dhi 


dpi 


< aa 


dhn 


is transferred through link q. If 


, then some energy 


dhi 


dpi 


> tta 


dhr, 


dprr 


then some energy must be called back, as explained in 
Section IIII-CI 

4) Go back to step 1, or terminate if sufficiently many 
iterations are performed. 


We describe our Algorithm in tabular form as Algorithm 2 
below. We note that our algorithm reduces to the one in ifT^ 
in the case of no energy harvesting or energy transfer. Next, 
we discuss the convergence and optimality properties of our 
algorithm. 


Algorithm 2 Algorithm to solve joint capacity and flow 
assignment problem for single time slot 

Initialize 

1 : Generate initial point 

Energy management step 

2: for n = 1 : iV do > All nodes 

Eind argmax,go_j(„) |^, perform (I?!]) as long as pi > 
— 1) is still satisfied 

3: end for 

Data routing step 

4: for n = 1 : do > All Nodes 

Eind path iF* that maximizes and ty* that minimizes 
X; where d G Odin) 

n,d 

5: for I € F* do = tf — 

6 : end for 

7: for / G (/* do tf^^ = tf +^t as long as pi > — 

1) is still satisfied 

8 : end for 

9: end for 


Energy routing step 


10: for q = 1 : Q do t> All energy links 

11: Set {i,j) ^ (origin,destination) of energy link q 

12: Set Ai = 


dhi 


dpi 


and A, = 

J dpj 


13: Use steps 6 : 10 of Algorithm 1 

14: end for 


15: Repeat until convergence 


tf^^ = tf and from (f7^ . this is only possible when ^ ^ 

is constant over all paths, which is equivalent to ItTI) . Using 
a similar argument we conclude that energy transfers satisfy 
(l72l i. This means that (t*,p*,y*) satisfies Lemma|9] ■ 

Now, we remark that even though we cannot claim global 
optimality of the solution, we have the following Pareto- 
optimality condition. 


B. Convergence and Optimality Properties of the Proposed 
Algorithm 

Every iteration of the algorithm decreases the objective 
function and the iterations are bounded. Using the fact that 
real monotone bounded sequences converge, we conclude that 
the algorithm converges. Assume (t*, p*, y*) is a convergence 
point of the algorithm. Next, we show that this point satisfies 
the KKT optimality conditions stated in Lemma |9l 

Lemma 10 (t*,p*,y*) satisfies the conditions stated in 
Lemma [9] 

Proof: When the algorithm converges, we must have pf^^ = 
pf. Erom (IT^ . this is only possible when ^ is constant for 
I G Odin) which is equivalent to (iTOl i. Similarly, we must have 


Remark 1 Assume that (t*, p* , y* ) satisfies the conditions 
stated in Lemma [9] then the vector of link delays is Pareto- 
optimal, i.e., there does not exist another pair of feasible 
allocations (t,p,y) such that 

hiipiji) <hiipf,tf), V/ (75) 

with at least one inequality being strict. 

This remark means that at the Pareto-optimal point, the 
average delay cannot be strictly reduced on one link without it 
being increased on another. The proof of this remark follows 
similar lines as the proof in IfT^ Thm. 4] and is omitted here 
for brevity. We note that, in particular, any local optimal point 
is Pareto-optimal due to the fact that local optimal points 
satisfy KKT conditions in Lemma |9l 
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Fig. 2. Network topology 1. 

VI. Numerical Results 

In this section, we give simple numerical results to illus¬ 
trate the resulting optimal policies. We study three network 
topologies shown in Figs. |2] [3 and |4] For all examples, we 
assume ai = 0.1 units \/l. The slot length is of 1 unit for 
convenience, so that we use power and energy; rate and data 
interchangeably. 

A. Network Topology 1 

We first consider the network topology in Fig. |2] with 
one source, one destination and three relays in between. 
The data and energy links are shown and labeled as 
in Fig. 12] where Z^s represent data links and y^s rep¬ 
resent energy links. The fixed data flows are t = 

= [2,1,0.5,0.125,2.125,0.375,0.5] units. We 
consider two time slots. The energy arrival vector is E = 
[(£;n,Ei2 ),..., (E4i,E42)] = [(15,10), (8 ,6), (5,9), (1,6)] 
units and energy transfer efficiencies are a = [ai , a 2 , 0 : 3 ] = 
[0.6,0.5,0.5]. 

The optimal energy transfer vector is found as 
y = [(2/11,2/12), (2/21,2/22), (2/31,2/32)] = [(0,3.75), 

(3.93,9.52), (2.35,9.81)] units and power allocation vector 
after energy transfer is p = [(pii,pi 2 ),..., (p 7 i,F 72 )] = 
[(7.5,7.5), (3.13,3.13), (0.62,1), (0.13,0.22), (9.17,11^ 
(0.45,0.74), (0.48,0.73)] units. Lemmas [6] and [7] can be 
verified numerically: h[{pii) equalizes for different outgoing 
links of the same node, for example, on links li and I 2 
(Lemma O; and where some energy is transferred, h[{pii) is 
proportional to the energy transfer efficiency of that energy 
transfer link, for example, / 12 ( 7 * 22 )/^ 3 ( 7 * 32 ) = ol\ (Lemma|2]l. 
Lemma [8] can also be verified numerically: after the energy 
transfers, the sum powers of the links are the optimal single¬ 
link powers. For example, node 1 has harvested (15,10) 
energies and transferred (0,3.75) of them. Equivalently 
node 1 has harvested (15, 6.25) and the single-link optimal 
powers for these harvests are (10.625,10.625) which is 
(7*11 + 7*21,7*12 + Vn)- It is interesting to note that node 



Fig. 3. Network topology 2. 

4 has transferred more energy than it initially had, which 
means that most of the transferred energy has been routed 
from other nodes. This is due to the high data flow on link 
Z 5 which leads to a higher energy demand at node 2 . 

B. Network Topology 2 

We next consider the star topology in Fig. [3 where 
five sources are communicating with one destination simi¬ 
lar to a multiple access scenario. The data flows are t = 
[0.5,2,0.5,0.5, 2] units. We consider a single time slot. The 
energy arrivals to all the nodes are the same, i.e., = 15 

units, Vn. The wireless energy transfer efficiencies are aq = 
0.5, Vy. 

The optimal energy transfer vector is found as y = 
[11.92,0,9.66,16.29,0] units and the power vector after en¬ 
ergy transfer is p = [3.07, 20.96, 5.33, 3.53, 23.15] units. This 
system is symmetric in terms of energy arrivals, channel noises 
and energy transfer efficiencies, and furthermore ti =t^ = 
and t 2 = fs. In this scenario, one might expect pi = P 3 = Pi 
and p 2 = P 5 - However, in the optimal solution p^ > P 2 - The 
reason for this asymmetry is as follows. Due to the high data 
loads on links I 2 and there is no incentive for these nodes to 
share their energy. Then, in the optimal solution, 2/2 = 2/5 = 0 
and nodes 2 and 5 act as energy sink nodes where energy 
is collected and not sent out. We see that node 5 has two 
nodes transferring energy to it while node 2 has only one 
node transferring energy. Then, p^ > P 2 - 

C. Network Topology 3 

In this last numerical example, we demonstrate the joint 
optimization of flow allocation and capacity assignment. We 
consider the diamond network topology shown in Fig. |4]where 
one source is communicating with one destination with two 
relays in between. The only exogenous data arrival to the 
network occurs at node 1 with the amount t = 2 units. The 
energy arrivals are [Ei, E 2 , E 3 ] = [2,0.5,1.5]. Energy transfer 
efficiencies are given as ai = 0:2 = 0.8. In this topology, there 
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are six unknowns to be determined, i.e., pi,p 2 , ^ 2 , 2 /i 7 2 / 2 - 
By exhaustively searching over these parameters, we can 
obtain the minimum achievable delay region as shown in 
Fig. EJtop). In the diamond network, there are two paths of 
data flow. One is the top path which includes links li and I 3 
and the other is the bottom path which includes links I 2 and 
I 4 . In Fig. Httop), we plot the delay on bottom path versus 
the delay on top path. Any delay which is to the interior 
of this curve is achievable whereas other delays are not. All 
points on this boundary are Pareto-optimal points. We observe 
that energy cooperation enhances the achievable delay region. 
In Fig. |3bottom), we demonstrate the convergence of our 
algorithm to a Pareto-optimal point. We start our algorithm 
from two different initial points and observe that they converge 
to a point which is on the boundary of the achievable delay 
region, demonstrating Remark [T] 




Fig. 5. (top) Achievable delay regions with and without energy cooperation, 
(bottom) Convergence of our algorithm. 


Appendix A 
Derivation of (1321) 


VII. Conclusion 

We considered the energy management and energy routing 
problems for delay minimization in energy harvesting net¬ 
works with energy cooperation. In this network, there are data 
links where data flows and energy links where energy flows. 
We determined the jointly optimal data and energy flows in the 
network and the energy distribution over outgoing data links at 
all nodes. We established necessary conditions for the solution, 
and proposed an iterative algorithm that updates powers, data 
routing and energy routing sequentially and converges to a 
Pareto-optimal operating point. In the special case of fixed data 
flows and no energy cooperation, we showed that each link 
should allocate more power to links with more noise and/or 
more data flow. In the case with multiple energy harvests, 
and no energy cooperation, we showed that the optimal sum 
powers on the outgoing data links of each node at every 
slot must be equal to the optimal single-link transmit powers. 
Our numerical results indicate that when data flows are fixed, 
energy is routed from nodes with lower data loads to nodes 
with higher data loads; while in the more general problem, 
where data flows are optimized also, allocation of data and 
energy flows are performed in a balanced fashion. 


Starting from OTl i. we have 






(76) 


We let r; = i log — ti, then 1 + ^ = With 

these definitions, we rewrite ( l76b : 


Or equivalently. 


A 


-2{ri+ti) 

2ai ' 


ne 



(77) 


(78) 


From here, r; = W{zi) where zi = and W{-) is 

the Lambert W function defined as the inverse function of 
w —>• we^ IIJTI . From the definition of rj, 

f log -f = r; = 1V(Z;) (79) 

and 

Pi = - l) (80) 

which is (I 32 I 1 . 
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Appendix B 
Derivation of (l34b 

From (l32l i. we have 

- l) 


(81) 


with zi = ^ 




. Our aim is to find To this end, define 


2\<7i 


d<7\ 


V / UZi OCi 


vi = — 1. Now we have, 

dpi 
dai 

The first partial derivative on the right hand side of 

— = p2(W(z,)+‘l)9_ ^(^0 

dzi 


zi{l + W[zi)) 

where we have used (l3Tl l. The second partial derivative in 


(82) 

is, 

(83) 


IS, 


dzi 

dai 




2Act/ ai 

Using dS^ and (l84l i in dS^ . we have 


ifi 

2 ai 


dpi 

dai 


= g2t, 


1 + W[zi) 


- 1 


(84) 


(85) 


which is (Ell. 


Appendix C 
Derivation of (1351) 

Starting from d32l i. we have 

pi+ai=aie^^'^^^d+U) (gg) 


with zi = Taking logarithm 

of d86] ), and differentiating both sides with respect to ti, we 
have 


1 1 dpi _ dWjzi) dzi ^ 

2 ai+pi dti dzi dti 

is evaluated from d^ and ^ is 


dzi _ 
dti 

Using (O and 

dpi 



— 2ti 


■pi) 


2 \ai 

in (Eli, we obtain 

W{zi) 


= zi 


— - 1 
2 ti 


= aie 


1 + W{zi) 

2 {W{zi)+ti) W{zi) -I 


— -1 
2 ti 

2 ti 


+ 1 


f/(l + W{zi)) 


which is 


(87) 


( 88 ) 


(89) 

(90) 


Appendix D 
Proof of Lemma[8] 

Assume that sum powers at each slot Si = bn is given 
for each i. Consider the inner optimization in dSSl l for a fixed 
slot, say slot i. For convenience, we drop the slot index i, and 


denote Si by s, and bn by bn We define a function x{s) as the 
minimization over bi for fixed s as follows: 



s.t. ^ 6; = s, 6/> 0, VZ (91) 

i 


which is the inner optimization in dSSl l for fixed i, and is also 
equivalent to: 


x{s) = min 

bi 


s.t. 



J2bi<s, bi>0, VZ (92) 

i 


Now, we claim that a;(s) is non-increasing and convex in 
s. Since increasing s can only expand the feasible set, x{s) is 
non-increasing in s. To prove the convexity: Let si, S 2 G 
Let 0 < A < 1 and A = 1 — A. Let bi be the solution of 
the problem with si, and b 2 be the solution of the problem 
with S 2 - Note that bi and b 2 exist and are unique due to 
convexity. The vector Abi + Ab 2 is feasible for the problem 
with Asi -I- As 2 since the constraints are linear. Then, 


a:(Asi -I- As 2 ) < ^ 

i 


_ ti _ 

1 log _ ti 

_AZ^_ 

1 log + -ti 


Xti 


i log + !^'j-ti 
= Xx(si) + Xx{s 2 ) 


(93) 


(94) 

(95) 


where d93b follows because the minimum value of the problem 
can be no larger than the objective value of any feasible point, 
d94b follows from the convexity of and (|95]) follows 

from the fact that bi solves the problem with si and b 2 solves 
the problem with S 2 - Now, the optimization problem in dSSl) 
can be written as: 

T 

min N x{si) 

Si 

i=l 

k k 

s.t. Gn Vi, Vfc (96) 

i=l i=l 

The problem in d96] l is in the same form as the problems in 
ifTrl eqn. (2)], ifTSl eqns. (6)-(8)] and ifTTl eqn. (15)] and is 
equivalent to the problem in lfT3l eqn. (3)], where a concave 
non-decreasing function of powers is maximized subject to 
energy harvesting constraints. In addition, m, ESI, CD have 
additional finite battery constraints which we do not have 
here. References ini, HI showed that the solution to this 
problem is invariant to the specific form of the function as 
long as it is convex (in minimization problems) or concave 
(in maximization problems). We follow the proof in HU 
Appendix B] and conclude that s, the optimal solution of ( 1961 ). 
is given by the single-link optimal transmit powers. 
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Appendix E 
Proof of Lemma[9] 

We show that the conditions in (l70b-(l72b are equivalent 
to (I62]i-(l64li therefore proving the necessity statement of the 
lemma. 

1) Writing (l62]) for node n and the data links I € Od{n) 
connected to it 


2ai 






< A„ 


Pi 

(97) 


Now, we claim that when pi > 0, Pi = 0. Assume pi > 0 
and Pi > 0. From dhST l. this means that pi = ai{e^*‘ — 1) 
and the delay at link I becomes ^ which is unbounded for 
ti > 0. Then, we must have ti = 0, but this means pi = 0, as 
otherwise power has been consumed on a link with zero flow. 
This is a contradiction to pi > 0. Thus, Pi = 0 when pi > 0 
and (ITOI i is satisfied with equality. 

2) We choose any origin destination pair (n, d) and identify 
a path starting from node n and ending at destination node d, 
and in which all link powers and therefore flows are strictly 
positive. We denote this path by J-n,d- We write the conditions 
(l63t on links on this path and sum them to get 







= E “ 2/3icrie^*' + 7z (98) 

Ti,d 

— ^ ^ l^n{l) (99) 

= Vd-Vn (100) 

= -Vn (101) 


where (l9^ follows from Pi = ji = 0 since pi > 0, ti > 0, 
(1 1001) follows from telescoping the sum J2ieJ^ ^ i^n(i) ~ 
and (IIOII) follows from setting Vd = 0 since it is a destination 
node and there are no flow conservation constraints at that 
node. We let = —v^ and get dZB. 

3) For energy link q between nodes n and m, k{q) = n and 
z{q) = m in d^ . From d^ . we have A„ = aqXm + dq > 
aqXm since 9q > 0. Then, 


ti 

2(71 







— 0(.q 


2,(Jk 






-1 


( 102 ) 

(103) 

(104) 


where dl021 i and (11041 ) are from using part 1 of Femma |9] for 
node n and m, respectively. Equality is achieved when pq > 0, 
since in this case 9q = 0 from (IFtI ). 
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